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Abstract Different initial and boundary value problems for the equation of vibra- 
tions of rods (also called Fresnel equation) are solved by exploiting the connection 
with Brownian motion and the heat equation. The analysis of the fractional version 
(of order v) of the Fresnel equation is also performed and, in detail, some specific 
cases, like v = 1/2, 1/3, 2/3, are analyzed. By means of the fundamental solution 
of the Fresnel equation, a pseudo-process F(t), t > with real sign- varying density 
is constructed and some of its properties examined. The equation of vibrations 
of plates is considered and the case of circular vibrating disks Cr is investigated 
by applying the methods of planar orthogonally reflecting Brownian motion within 
Cr. The composition of F with reflecting Brownian motion B yields the law of 
biquadratic heat equation while the composition of F with the first passage time T t 
of B produces a genuine probability law strictly connected with the Cauchy process. 

Keywords: Schrodinger equation, higher-order heat equations, Brownian motion, Mittag-Lcfflcr 
function, Fresnel function, Wright function, inversion radius, clastic Brownian motions, fractional 
diffusions, vibrations of plates. 
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1 Introduction 

One of the most important equations of mathematical physics is that of the vibrations of rods (or 
equivalently of beams and shafts) which can be written as 



a u 
W 2 



x £ R, t > 0. 



(1.1) 



The constant n appearing in (1.1) is related to the physical structure of the vibrating rod and will be 



considered below as equal to k = 1/2 for reasons which will appear clearly further in the text. For the 
derivation of (fOJ consult [5J pages 112 - 114] or the classical book by Courant and Hilbert [5J pages 
244 - 246] where, also the derivation of the equation of the vibrations of plates (pages 250 - 252) 



dt 2 



dx 2 



dy 2 



(1.2) 



is presented (p is a physical constant assumed below equal to p = 1/2 2 ) 
written as 

d i d 2 \ / d id 2 
dt + 2~lh?) (dt " 28X 2 " = !1 



The equation (1.1) can be 
(1.3) 



and this shows a strict connection with the Schrodinger equation and therefore with the heat equation. 
In effect 

du i d 2 u 

dt =± 2dx 2 (L4) 
can be easily reduced by means of a time change t' = ±it to the classical homogeneous heat equation. 



This connection entails that initial-boundary value problems concerning equation (1.1) have solutions 



which can be constructed by means of related solutions of problems concerning the heat equation and 
therefore are related to the standard Brownian motion B. In particular, the fundamental solution of 



( 1.1 ) can be obtained by means of the rule 



u(x, t)dx =- 



Pr{B(s) £ dx} 



dx 

"~2~ 



2izt 



'2lTt 



Pr{B{s) £ dx] 
dx 



- [p{dx, t} + p{dx, —it}] 



V2irt 



cos 



.r 
2t 



(1.5) 



The measure p appearing in (1.5) must be understood in the sense that 



p{dx,±it} = Pr{B(s) £ dx} 



s—±it 



The idea underlying ( 1.5 ) is that at time t — a Brownian motion is started off from x = either with 



increasing or decreasing imaginary time. The choice of time direction is made once only and cannot 
be changed. This process can be regarded as the limit of a random walk with symmetric real-valued 



steps separated by imaginary time intervals. Each term in (1.5) is therefore connected with a sort of 



Brownian motion where time takes imaginary values. This approach permits us to exploit the panoplie 
of mathematical results and instruments of Brownian motion to obtain solutions of all possible types of 
problems for vibrations of infinite, semi infinite and finite rods. Furthermore, the finite, signed density 



of the measure (1.5) suggests the construction of pseudo-processes (called Fresnel pseudo-process and 
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denoted by F(t), t > 0) whose finite-dimensional distributions can be constructed according to the 
rule 



1t(X] 



n 1 
, i i . . . , x n , £ n ) J^J — — 



Pr{B(si) e dxi, B{s„) e dx n } 



j = 1, 2, ... , n. The construction of Wiener-type measures for higher-order heat-type equations 



du 
St 



±- 



S n tt 



(1.6) 



(1.7) 



has been realized by different authors since the beginning of the Sixties (see Krylov |12] . Ladokhin 
[16]). Much work on pseudo-processes connected with (1.7) has recently been done by Lachal (15) 113]. 
Cammarota and Lachal [15] . where by different means and techniques, various types of functionals of 
these pseudo-processes have been investigated. 

The fractional version of the equation of vibrations of rods 



d 2v u 



i cru 

¥ Sx 1 ' 



is also examined. The Fourier transform of the solution of (1. 
conditions, takes the form 



< v < 1 (1.8) 
equipped with the necessary initial 



U 2u (J3,t) 



+ OG 



A/3x 



U2v{x, t)dx 



E„ 



./3 2 t L 



,P 2 F 



= E 2 



$H 2v 
2 2 



where 



E Vil (x) 

is the Mittag-LefHer function. We are able to invert the Fourier transform above and obtain an explicit 



^ T(vk + 1) 

fc=0 v ' 



solution of the fractional equation (|1 
1 



U 2 v{x,t) 



OC 



i 



m=0 



V2 



in the form 



cos 



1 



1 



r 



i 



(1.9) 



From (1.9) for v = 1 we obtain the fundamental solution (1.5) while for v = 1/2 we get one alternative 
expression for the solution of the biquadratic heat equation. This permits us to give a sketch of the 
function ui(x,t), x € K., t > which looks like the normal bell-shaped density in the neighborhood of 
the origin (see figure [6]). In great detail we examine the function U4/3 (y — 2/3) and prove that it is 
possible to reduce it to the superposition of Airy functions. Also, the case v = 1/3 is analyzed and we 
prove that it can be expressed in terms of (1.5 1 and the Airy function. 



The analysis of Section 3 shows that the profile of the vibrating rod has a peak near the origin 
where the initial disturbance is originated and symmetric damping waves which rapidely decrease in 
size. In our opinion this accords with what happens with real vibrations where friction dissipates the 
original energy and therefore fractional equations of vibrations of rods better describe the phenomenon 
under investigation, see figure [H] 



Section 5 is devoted to the multidimensional version of the equation ( 1.1 ) which governs the vibra- 
tions of plates. Its general form is 



dt 2 




(1.10) 
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and possesses the fundamental solution of the form 




u(xx,...,x d ,t) = . cos [ ~ d J ) ■ ( L11 ) 

In the plane we study the vibrations of a circular plate Cr with the Neumann boundary condition on 
the edge 8Cr. We solve this problem by applying an approach based on the inversion with respect to 
the circle which parallels the construction of the reflecting planar Brownian motion inside a disk. This 
is useful to describe the oscillations of thin vibrating structures, started off by a concentrated central 
initial disturbance. 

Finally, we show that the composition of the Fresnel pseudo-process F(t), t > with reflecting 
Brownian motion |-B(t)| produces a subordinated process F(\B(t)\) whose one-dimensional law coin- 
cides with the fundamental solution to 

du 1 d 4 u 

dt ~ 2 3 ftz 4 ' fl 12) 

u(x,0) = 6(x). { ' ' 

Instead, the composition of F with T t = inf{s > : B(s) = t} yields the following interesting genuine 
probability law 

Pr{F{T t ) € dx}/dx = -^ t 4 + X 4 , x£R,t>0 (1.13) 

7TV 2 t + X 

which solves the fourth-order equation 

d 4 u d 4 u „ 

^ + ^=° (L14) 



The distribution (1.13) has two maxima which move in opposite directions as time passes and display 
a structure similar to the solutions of fractional diffusion equations for a degree of fractionality varying 
in 1 < v < 2 (see [S], [Hj). Finally the successive compositions of Fresnel pseudo-processes and the 
law of 

F n (t)=F 1 (\F 2 {\...F n+1 {t)...\)\), i>0 (1.15) 
is shown to be governed by the higher-order equation 



w ( x ,t) = -2-^ -D^^i). (1.16) 

This is clearly related to what happens with the n-th order iterated Brownian motion as shown in 
Orsingher and Beghin [21) . or also by Baeumer et al. PQ for fractional equations. For higher-order 
equations, see also D'Ovidio and Orsingher [7]. 

1.1 Notations 

• u(x, i) = -^=| cos — j J : fundamental solution of the one-dimensional Fresnel equation, 
. U(fi, t) = f+™ e^ x u(x,t)dx, 

• u a is the solution to the Fresnel equation with absorbing conditions, 

• u r is the solution to the Fresnel equation with reflecting conditions, 

• u el is the solution to the Fresnel equation with elastic boundary conditions, 

• p el is density of the elastic Brownian motion, 
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• q is the fundamental solution to the biquadratic equation. 

• p re f (r, t) is the probability density of the reflecting Brownian motion in a disk Cr (q re ^ is its 
kernel) , 

• p re f is the law of a vibrating plate of radius R with Neumann condition on the border OCr (q re ^ 
is its kernel), 



• U2 U is the fundamental solution of the fractional Fresnel equation (1.8), 



v v is the fundamental solution to the fractional diffusion equation (3.17) 



The general unknown functions of the equations appearing in the paper (also the fractional equations) 
are indicated by u while the solutions of the boundary Cauchy problems are denoted by u a , u el , etc.. 
The solutions of the initial-value problems for fractional Fresnel equations are denoted by Ui v and for 
the fractional diffusion equation by v v . 

• u(xi, . . . , Xd, t) = ^y^y cos (y2j—i ^| — d^j : fundamental solution of the n-dimensional Fres- 
nel equation, 

. Ufa, ...,p d ,t)= /+~ dxi... f+™ dx d e^U ^u(x u ...,x d , t), 

(27r)-"/ 2 (^(Xi-Xi-l) 2 TT 

u(xi, t\\ . . . x n , t n ) — === - cos > — r n — 

K ' n"=i vtj - tj-i ypiWi-ti-i) 4 

is the n-dimensional measure density of the Fresnel pseudo-process F(t), t > 0, 

• B{t), t > is a standard Brownian motion, 

• F(t), t > is the Fresnel pseudo- process, 

• T t = inf{s > : B(s) = t}. 

2 Vibrations of rods and Brownian motion 

The equation of vibrations of a thin, flexible rod has the following form 

d 2 u 1 d 4 u m 

Vibrations of rods differ substantially from vibrations of strings (and membranes in dimension 2) 
because strings propagate with infinite velocity as heat does. Furthermore, the heat diffusion and the 
vibrations of rods have strict connections which inspire much of the work of this section. 



Since (2.1) can be written as 

d i d 2 \ ( d i d 2 \ 
we can decouple the original equation into two Schrodinger equations (each of which can be reduced 



to heat equations by a suitable change of variable, t' — ±it). The fundamental solution to (2.1) can 
be obtained by suitably combining the solutions to 
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The fundamental solutions to (2.3) and (2.4) can be written as 




e ' 4 
\f2wi 


COS 


x 2 
21. 


+ isin 


X' 2 
21 


V2TTt 


COS 


2 

X 

21 


— i sin 


x 2 
21 



and thus 



u(x,t) =- vi(x, t) + v 2 (x, t) 



1 



'Ant 



2f 



2t 



'2irt 



We call (2.6) Fresnel function and can easily check that 

/ u(x, t)dx = 
Ja. 

because 



cos x 2 dx 



sin x 2 dx 



1 TT 

2V2- 



The function (2.6) is the solution to the Cauchy problem 



d 2 u _ _J_cfu 
dt 2 ~ 2 2 dx 4 

u(x, 0) = 8{x) 
u t (x,0)=0 



and its Fourier transform reads 



1709, f) 



3 tl3x u(x,t)dx = COS 



f3 2 t 



This can be proved by writing 

f+OO 



1 



e i0x cos 



x 
2t 



— | dx = -e 
4 / 2 



/.+oo ,-<£±±ir 



+ -e 4+ 2 



+oo 



V2nt 



-dx 



12-nt 



-dx 



x 
2f 



2^2 p 2 t f°° 2 , p 2 t 
= — =• cos — — / cos y dy — cos — — . 



Result (2.9) can be also obtained by solving the Fourier transform of problem (2.8), that is 



a 2 u _ 

dt 2 ~ 

t/(/?,0) = l, 

W,o) = o 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



The structure of the function (2.6 1 is illustrated in figure [Tj and it represents the profile of the vibrating 
rod at different times. It shows that the deformation instantaneously propagates along the x axis and 
the oscillations decay as time passes. The measure density ( |2.6[ ) is a sign- varying function and plays an 
important role in the theory of optical diffraction. By studying diffusion processes with an alternating 
time direction (forward or backwards) an equation similar to (2.2) was obtained and its solutions 
investigated in Orsingher [2"0] . 
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Figure 1: The Fresnel function (1.5| represents the propagation of vibrations along an infinite rod at 
time t = 1,20,40,60. 



Remark 2.1. The area between two successive roots a k , ctk+i of the Fresnel function (1.5) tends to 
zero as k — > oo. This can be shown by observing that 



1 



/2vrf 



COS 



It 



dx 



< 



1 



X IX 

. . cos — 

a k \2rrj ak t \2t 



dx 



< 



1 



~a k V 2ir 



sin 



2f 



— sm 



fc+i 



2t 



where a k — \/2TTyJ?>/2 + k. 

By exploiting the ideas outlined above it is possible to study the vibrations of semiinifinite rods 
with different types of boundary conditions. We start by taking into account the vibrations with an 
absorbing condition at the free end point x = 0. The rod is assumed to be clamped at infinity. We 
arrive now at our first Theorem. 

Theorem 2.1. The solution to the boundary-value problem 



SsKM) = z>0, t>0 
u(x, 0) = S(x — y) 
u[0,t) = 
ut(0,t) = 



is given by 



u a (x,t;y,0) 



2irt 



(x — y) 2 7T 
2t 1 



2t 4 



(2.11) 



(2.12) 
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Proof. Our scheme consists in considering two absorbing Brownian motions, one with imaginary in- 
creasing time and one with imaginary decreasing time. Thus, we can write 



u a (x, t; y, 0) = - [p(x, it; y, 0) + p(x, -it; y, 0)] 



(2.13) 



where p(x, s; y, 0) is the transition function of a standard Brownian motion. Each component of (2.131 
solves the boundary- value problem for the Schrodinger equation 



du 
dt 



±i 



dx 2 



From (2.131 some calculations lead to (2.12). 



□ 



We can observe that (2.12) solves the equation of vibrations of rods and the boundary conditions 
can be checked by setting x — 0. In the same spirit we can solve much more complicated boundary- 
value problems such as those involving the elastic condition at the free end point of the vibrating rod. 
Elastic Brownian motion was studied by Ito-McKean jTSj . The fractional extension of elastic Brownian 
motion was investigated by Beghin and Orsingher [2]. Our result is in the next Theorem. 

Theorem 2.2. The solution to the elastic boundary-value problem 



Qf_2 (^5 ^) — 2 2 dx A X > 0, ^ > 

u(x, 0) = S(x — y) 
u t (x,0) = 

u x — au\ x= Q = 0, a > 



(2.14) 



is given by 



u e \x,t;y,0) 



1 

^2kI 

- 2e a( - x+v ^ 



(x - y) 2 _ 7T 
2t 4 

+oo we -aw 
V2-Kt 3 



1 



'2-Kt 
2ae a{x+y) 
V2~Tri 



x+y 

(x-y) 2 

2t 



— cos 

~2t 
cos 
w 2 



(x + y) 2 
2t 
3tt 

" T 



dw 
(x + y) 2 



x+y 



2t 



2t 
dw. 



(2.15) 



(2.16) 



Proof. The solution to (2.14) can be written as 



u el (x,t;y,0) = X - [p el (x,it;y,0)+p el (x,-it;y,0)] 



(2.17) 



where p el (x, s; y, 0) is the solution to the heat equation with the elastic condition at x = which can 
be written as 



p el (x,s;y,0) 



y/2TTt 



'2-Kt 



It 

J2wt 

it 

J2wt 



2 £ a(x+y) 



2ae a(x+y) 



we 



x+y 



e 2t 
- aw ^=dw 



e 2t 

--dw 



x+y 



'2-Kt 



(2.18) 

(2.19) 



By considering (2.18) and (2.19) we arrive at the explicit solution of problem (2.14) which coincides 
with (|2.15[) (or alternatively with (|2.16|)). □ 
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Remark 2.2. For 



oo from (2.181 we extract the absorbing solution (2.11). For a = we obtain 



the reflecting solution, that is the reflecting solution to the equation of vibrations of rods with the 
Neumann condition at x = 0. We can extract from (2.18) the second expression (2.19) by means of 
an integration by parts. In order to check that the solutions (2.15) and (2.16) satisfy the equation of 
vibrations of rods it is convenient to rewrite the third term of ( |2.16 ) as 



2a 



: COS 



((x + y) + w) 2 
2t 



— - dw 



Remark 2.3. We note that the boundary conditions 

tt(0,f) = 

and 



dx 2 



(M) 



(2.20) 



lead to the solution (2.12). The condition (2.20) means that the bending moment at the free end point 
x = vanishes. The "reflecting conditions" 



du 



and 

both lead to the solution 

u r (x,t;y,Q) = 



Ox 
dx 3 



12-nt 



cos 



(x,t) 

x,t) 

(x-y) 2 
2t 







x=0 



COS 



2t 



Remark 2.4. We now evaluate the "survival measures" of the elastic, reflecting and absorbing kernels 
obtained so far. We restrict ourselves to the elastic case since the other two follow without effort. 



u el (x, t; y, 0)dx =- 



1 



'2-Kt 

- 2a 
1 

y/2M 

- 2a 
1 

\/27rf 



cos 



e a(x+y ^dx 



2t 



x+y 



COS 
o 

w 
2t 



(x + y) 2 



2t 



dx 



dw 



cos 



w 
~2t 



j dw + 


f°° fw 2 




/ cos — 






Jy \2t 





dw 



fw 2 






-;)*jT 


\2l 





e <*(*+y) dx 



cos 



w 
~2f 



dw + 2e c 



w 
~2t 



dw 



(2.21) 



For a = we have the reflecting case and the integral above equals one while for a — ¥ oo the second 
term in (2.21 ) goes to zero and we retrieve the "survival measures" of the absorbing case. 

For rods of finite length we can get explicit solutions by resorting again to the decomposition of 
equation (2.1) on considering the solution of the corresponding Schrodinger equation. We have that 
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Theorem 2.3. For the rod of finite length L we have that the solution of 



2 2 a 

u(x, 0) = 8{x — y) 
u t (x,0) =0 

, U x \ x =0 — Ux\x=L = 



is given by 



+oo 



/c — — oo 

Proof. The solution can be written as 



(x-y + 2k,L) 2 7T 
2t 4 



(x + y + 2kL) 2 7r 
2t 4 



(2.22) 



(2.23) 



u r (x, t; y, 0) = - it; y, 0) + -it; y, 0)] 



where u(x, t\ y, 0) is the well-known solution of the Cauchy problem of the heat equation 



ffOM) = 5§J0M), 0<x<L, i>0 



u(x, 0) = 5(x — y) 
u t (x,0) = 

U x \ x -o = U x \x=L = 0. 



(2.24) 



□ 



3 Fractional equations of vibrations of rods and fractional dif- 
fusions 

3.1 General results 

In this section we consider the fractional version of the equation of vibrations of rods 

d 2v u 1 d A u 

< i/ < 1, x6l, t>0 (3.1) 
u(x ) 0)=5(x), 0<i/<l (3.2) 



dt 2 " ' 2 2 dx i ' 
subject to the initial conditions 



and also 

wt(a;,0) = for 1/2 < v < 1 (3.3) 



The differential operator appearing in (3.1) must be understood in the sense of Dzhrbashyan-Caputo, 
that is (for information on this derivative consult Podlubny [2~4"j ) 

BTu t _ ±s 1 f* %£(x,s) 



(x, t) = — -. r / . . — -. ds, for m — 1 < v < m, m S N. 

dt vX ' r(m - v) J Q (t - s)» +1 - m ' 



For v = 1 equation (3.1 1 coincides with the equation of vibration of rods and beams while for v — 1/2 
it gives the biquadratic heat equation. The latter equation has been investigated by many researchers 
since the Sixties and even a stochastic calculus related to it has been worked out (see for example 
Krylov [T2], Hochberg (TU], Nikitin and Orsingher [25], Lachal [HUH]). 



Our first result concerns the Fourier transform of the solution to ( |3.l[ ) (with initial conditions (3.2 ) 
and Q). 
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Theorem 3.1. The Fourier transform 

/+oo 
e if3x u v {x,t)dx 
-OO 

of the solution to the Cauchy problem 



(3.4) 



§t^ = -2W> 0<Kl,xel,t>0 

u{x, 0) = S(x) 

u t (x,0)=0 



reads 



where 



.P 2 t v 



E2V.1 



2 2 



k>0 



T{vk+l) 



(3.5) 



(3.6) 



(3.7) 



is the Mittag-Leffler function. 



Proof. The Laplace transform of the solution to (|3.5 ) reads 

,v-l 



111 / e l0x u 2lJ {x,t)dxdt 



1 



^ + if3 2 /2 ^~ifi 2 /2 



/' 



2i/-l 



H 2v + /3 4 /2 2 ' 



The inverse Laplace transform of (3.8 1 yields (3.6 1. 



(3.8) 



□ 



For the evaluation of the inverse Fourier transform the first expression of (3.6) is particularly 
convenient. This is because each term in (3.6) represents the Fourier transform of the solution to the 
fractional diffusion equation 



f^ = A 2 §^, 0<i/<l,i6R,{>0 



u(x, 0) = S(x) 
ut{x,0) = 0. 



(3.9) 



The solution to (3.9) can be written down in terms of Wright functions 



fc>0 



k\T{ak + P)' 



a > -1, /? > 0. 



We now arrive at the next result, concerning the inverse Fourier transform of (3.6). 



Theorem 3.2. The solution to (3.9) is given by 



tt\/2F ^ n ml \ y/F 

v rn— v v 



COS 



m + 1 



771+1 



■kv r 



m + 1 



'2t v 

m— 



y/¥J r(i 



cos(^Tr) 

m-j-1 



(3.10) 
(3.11) 



/or a; S M, i > and < ^ < 1. 
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Proof. For the fractional diffusion equation (3.9) the solution reads 

v v (x,t) 

and has Fourier transform 



1 

2\Jt 



V 2 < x 2 



e lf3x v„(x,t)dx = E u<1 (-A 2 /?V) . 



(3.12) 



Therefore, in view of (|3.6|) we have that 
1 



U2„(x,t) 



l 



l 



=W_ K -W — 



2 |z| 



■ 1*- 1 



2 > A 2 



2^_ 



-V2^e 



1 OO 

= E 



'2 3 t 



m— 



e 4 * ( -v^ 



|x|e l < 



-|a;|e l J 



m=0 
oo 



FJ m!r(-^f + l-|) 



v £- — t \ 



COS (^7r) 



& J m!r(-^ + l-|) 



v m— 



COS 



m + 1 



m + 1 



7r^ r 



m + 1 



In the last step we applied the reflection formula 

r(z)r(i-z) 



SUITTZ 



Remark 3.1. We give now the Laplace transform of (3.10). Since 



p oo 

Jo V 2 2 



and 



/ 1/ 




'V 




1 m — 


i) r ( 






V 2 




,2 mH 





sin^!/ 



we extract from (3.10) that 



f°° 1 °° 1 / \ 



H 2 m+ 2 1 C OS 



m + 1 



2^2 m! V 



/2 



I//2-1 



2V2 



M "72-i e -MM" /2 



2^ 



(3.13) 



□ 



12 



^«V2-i e -klft" /3 

: 7^ 



cos (W/i"/ 2 - J 



(3.14) 



We now take the Fourier transform of (3.14) 

r-l 



A* 2 



2v^2 

r-l 



i/3a= e — |a=lM" /2 e" 4 + e -if 



'V' 2 



e i/3a; e -|x| At , ' /2 e 



-j-OO 



,3 2 + 2^e l 2 /?2 + 2^e"' 
i(^ 2 - 2i/z I/ ) - i(/? 2 + 2^") 
2 V + /5 4 



A 4 + 2 2 



and this confirms result (|3.8|). For a further check we evaluate the inverse Fourier transform of (|3.8|). 

,+oo e - lf}x d p 



2i/-l r+oo -i/3x 



2n 



„2l/ _(_ /3 2 27T 



2 2 7T 



-OO /i^ ~T~ 



+oo —i@x 



-i(5x 



oo [if 2 

(2e-«fM") 1 / 2 2 



2 J 



d/3 



4-oo 



2V2 



(2e-^^) + /3 2 e'f (2e- l f^") 1 /2 
(2e 1 ^") + /3 2 e-*f (2e i t^ 1 ") 1 /2 

-|a;|(2e*f- M l ') 1 / 2 + i f - |s| (2e~ 4 5 ft" ) 1 ' 2 -i f 



/' 



2V2 

%-l e -\x\n" /2 ^ p 

71 2 



(2e 4 f^) 1 /2 (2e- i f^) 1 /2 



H(i-i)M" /2 -if + e -kl(i+«)M" /2 +*7 



e i|x|M" /2 -iJ _|_ e -i|x|/i" /3 +if 



V2 



COS 



(w^ /2 - \ 



Theorem 3.3. The solution to the fractional equation of vibrations of rods 

d 2 "" — A 2 d^u T rD /- •> f) 

2 2 C?X 4 ' ' 

u(x,0) = <5(x), < ^ < 1 
Ut(x,0) = 0, 1/2<i/<1 



coincides with 



u 2v {x,t) = 



1 



: COS 



V2^ V 2s 4 

where v 2v is the folded solution of the fractional diffusion equation 



v 2 u(s,t) ds 



subject to the initial conditions 



d 2v u X 2^ 2u 
~IW = dx 2 



u(x,0) = 6(x), 0<v<l, 
u t (x,0)=0, 1/2 <//<!. 



(3.15) 

(3.16) 

(3.17) 
(3.18) 
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Proof. The fractional derivative of order 2v of (3.16) yields 



d 2v u 2v 
dt 2v 



: COS 



x 7r\ d v 2v 



2s 4 

OO 



dt 2v 



(s, t) ds 



r°° i 

= [by pllf ] = A 2 / -— = cos 

JO \J1~KS 



2s 



ds 2 



(s, t) ds 



r°° d 2 


1 


fx 2 






. cos 




-;)] 


Jo M 


_V*KS 


v 2s 





v 2v (s,t)ds = [by (2.1)] = 



A 2 G> 4 ?x 2t , 
2 2 <9z 4 ' 



We note that 



/ e ifix u 2 y(x,t)dx = [by ([3~T6|] = / cos 
Jr Jo 



P 2 s 



V2v(s,t) ds 



I 

Jwt 



e 1 ^ v 2 y{s,t)ds = E 2l/ ,i 



A 2 /3 4 
2 2 



□ 



Remark 3.2. Theorem 3.3 is the counterpart of result (2.2) for fractional diffusions treated in [21] 
where the role of the Gaussian density is here played by the Fresnel function (1.5). In some cases 
this is extremely fruitful because it permits writing down the fundamental solution explicitely. For 

(3.19) 



v = 1/3, from (3.15) we obtain that 

f°° 1 fx 2 tt\ 

U2/s(x,t)=J - / ^cosl—--\v 2/3 (s,t)ds 



1 



/2ns 



7T\ 3 1 



2s 4 J A 



Ai 



ds. 



In the intermediate step of (3.19) we employ the result 

d 2 ' 3 u 



= A 2 |^, x£R,t>0 



g t a/3 — n dx 2 

u(x, 0) = 5(x) 



(3.20) 



whose solution is (see [2TJ formula (4.2)] 

"2/3 OM) 



3 1 



Ai 



iei,t>o. 



Formula (3.19) is an integral version of (3.10) for v — 1/3. The relationship (3.16) relates the solutions 
of the fractional equation of rods (3.15) with the fractional diffusion equation (3.17) of which a vast 
literature exists. 



3.2 Special cases 

We now examine some special cases. 



Remark 3.3. For v = 1, we are able to extract from (3.10) the fundamental solution (1.5). For n £ N, 
we consider in (3.10) the cases m = 2n + 1 (for which the series is equal to zero) and m = 2n for which 
we have that 



u 2 (x,t) = 



1 



1 



nV2l ^ ml 

v m— 



-V2 



t 1 / 2 



cos 



1 



1 



7T r 



1 
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Figure 2: The fundamental solution u 2 /z of the time-fractional Fresnel equation. 



n=0 

OO 



2n + 1 

cc» I 7T | Sill 



2n + 1 



7T r 



2n+l 



^2t^ Q (2n)! V ' j V 4 



^E 



'2-Kt 



E 



(-1)" f\x\ 



n=0 
oo 



n 7r\ r(2n) 



(2n)! V 2t 
I ^ (-1)" /^|a;| 2 V l 1 



/n 

— 7T 

V2 



4/ T(n) 



'2irt 



n=0 



„z -if 7T-M t i „— — »- 



1 °° 1 



2t J 2 

2 \ ™ 



2i 



2 \ ™ 



2t 



1 



g 1 2t + l 4 -j- g+* 2t * 4 



1 /re* 7T 

-7 cos — - 

/27rt l2t 4 



2V27ri 

Remark 3.4. Another important case is v = 1/2 and in this case we easily arrive at 



2tt 



e 22 cos (xy) dy 



2nV2t 



'—f^t-^MTi-irr 



k 1 

2 + 4 



(3.21) 



The integral (3.21 ) had been studied by Bernstein |4] and its connection with the fundamental solution 
of the biquadratic heat equation is well-known. By applying the multiplication formula of the Gamma 
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function 



r(»r [z + —)T[z+—)...T[z + 



m — 1 



(2K)^mi- mz T{mz). 



for to — 4 and z = fc/2, we have that 



r(| + |) _ n2i-™ 



r(2fe + i) fcir (| + 1) ' 



This leads to the alternative form 



ui(x 



1 V s 



2^^ Q V2^; fc!r(| + | 



(3.22) 



Particular attention is drawn to the case v = 2/3 for which we produce different alternative 
expressions of the solution in terms of Airy functions. 



Theorem 3.4. The solution to 



with 



d^u _ 1 d 4 u 

<% 4 /3 ~ ^22^4' 



ieR,t>0 



u(x, 0) = S(x) 
u t (x,0) = 



has the form 



where 



u 4/3 (x,t) 



2 3 / 2 s/3i 



e^Ai { V2 



\x\e * 



■e'^Ai y/2 l - 



Ai(x) — — I cos ( ax H I da 



"5 V3 



I ! -X^ \-lA -X 3 ' 



5 V3 



is £/ie ^4in/ function (the second formula in (3.26) holds for x > 0). 



Proof. For v = 2/3, formula (3.101 yields 



m—0 



m + 1 



m + 1 \ /m + 1 

7r r 



3 



We split the series (3.27) into the following three cases: 

TO = 3fc + 2, 

m = 3fc + 1, 
to = 3fc. 

For our convenience we write 

U4/3OM) = M + A a + A . 
For to = 2fc + 2 the series ^2 equals to zero because 

'to + 1 

sin " 



7T = sin((fc + 1)tt) = 



(3.23) 
(3.24) 

(3.25) 
(3.26) 



(3.27) 



(3.28) 



16 



for all integer values of k. For n = 3k + 1, formula (3.10) takes the form: 

3fc+l 



1 y> 1 

irVz^t ^ (3k + 1)! 



cos 



3k + 2 



7r sin 



3/j + 2 



7T r [k 



"^727^^ + 1)! v ^/ 



3fc+l 



3/s 



sin — 7r (-1) sin -ttT k 



(3.29) 



In view of the triplication formula (see [El page 14]) 



T(z)T z+- T z + - = 



2tt 



3 3z-l/2 



r(3 2 ) 



we have that 



1 2tt 



r(3fc + 2) 3 2 3 3fc - 1 /2r(/c + i)r(fc + i + i) 



(3.30) 



The series (3.29), in force of (3.30), becomes 

3fc+l 



_2 7 rsinf 1 



3k 

sin | — 7r 



3fc+l 



1 y 1 (rz \x\ \ 
72 7* k\ V Z¥t) r (k + § + 1) 2z 



1 1 

4 -y 3 3fc - 1 /2r(fc + i)r(fc + 1 + 1) 
i i r 



E 

fe=0 



e "4 



72- 



37* 

x I e l * 



3fc+l 



e' 4 



7^ 



37* 



3fc+l 



1 



72 



cc|e * 



3^i 

1/2 oo 



fc=0 



37* 



3/2~ 



2fe+i 



fc!r(fc + i + i) 
i 



fe=0 



37* 



3/2\ 2fc +3 



fc!r(fc + i + i) 
i 



V2 



1/2 



72 



37* 

1/2 

A/3 ( 2 ( V2 



A/3 2 ( 72 



|x|e 4 
37* 



fcir (fc + i + i) 

3/2\ 



x e 4 x 



37* 



2s/2^t 



-e" 4 



72 



x|e 4 



1/2 



A/3 2 ( 72 



72 



37* 



A/3 2 [ 72 



|x|e 4 
37T 

3/2 



3/2N 



With similar calculation we evaluate 



An = 



1 °° 1 

— T — 

m 3/7 rafc 



7r727i H (3fc)! V 7* 



fc=0 
gl/2 00 



cos 



37* 



3fc + 1 
4 

3fc + 1 



(3.31) 



tt sin 



3fc + l 



7t r 



3fc + 1 



tt r fc 
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31/2 



3 A: 



3fc + 1 



1 



E 

fc=0 



e" 4 



72 



x|e 4 
~3~W~ 



4 y 3 3fc + 1 /2r(fc- i + i) 



3 /11 - .» 2L. \ 3ft 

x e 4 



e ' 4 



V2i 



fc!r(fc-| + i) 



V2 



\x\e l 4 
3Vi 



EN 5 



fc=0 



|x|e 4 
3^ 



V2 



I 1 — t — \ l/ 2 00 

i e * 4 x 



3^i 



£ 2 V2 



1 



2 3 /2^ 



72 



xje* 4 



fc=0 

1/2 



3^t 



3/2* 



3/2* 



2fc-i 



fc!r(fc - l + i) 



2fc-i 



fcir (fc - § + 1) 



I-1/3 2 I V2 



|x|e l4 



3/2N 



e' 1 * ( V2 



3^ 



3^t 



By summing up A^+i and A^k we have 

1 



«4/a(a:,*) =7^/1 



V2 



\x\e 4 

"3?r 



1/2 



J_ 1/3 2 ( V2 



3^ 



3/2N 



-e' 1 * ( y/2 



-e" 1 * I \/2 



3^ 



23/2^ 



e*? \/2 



3^ 



e _i *Az [ 72 



(3.32) 



□ 



Remark 3.5. Result (3.251 shows that the solution to 

<9 4 / 3 u _ 1 d 4 u 

can be expressed as 



(3.33) 



"4/3 OM) =2 [ w 2/3(^,te l4 ) + u 2 / 3 (Me l4 ")] 



' ,1, / ' 



2 2 ^3t LAi VAiV3i/ A 2 VA 2 v / 3i 



with 



where 



Ai = 



A 2 = 



, \ 3 1 , 



e 4 
75 



2Av / 3t VAv 7 ^ 



(3.34) 



(3.35) 
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Figure 3: u 4 / 3 : the function (3.25) 



is the solution to 



dt 2 / 3 = dx 2 ' 



u(x, 0) = S(x) 



(3.36) 



(see formula (4.2) of 2lJ). Formula (3.34) shows that the solution to ( |3.33 1 can be viewed as the 
superposition of solutions of fractional diffusion equations of order 2/3 with imaginary time running 



in opposite directions. Result (3.34) can also be obtained by inverting the Fourier transform (3.6) for 
v = 2/3. 

Remark 3.6. By introducing the drift in the Schrbdinger equations appearing in the equation of 
vibrations of rods we get that 







d_ 
di 

dt 2 



dx 2 

1_ (&_ 
2 2 [dx 2 



d_ 

dx 



d_ 
di 



dx 2 



d_ 

dx 



dx 



The fundamental solution to (3.37) takes the form 

u(x, t) 



'2nt 



u. 



x 

Yt 



(3.37) 



(3.38) 



which slightly extends (2.12) 



4 Pseudo-Processes related to the equation of vibrations of 
rods 



We introduce a pseudo- process associated to the fundamental solution of equation (2.1 ) as follows. We 



can assume that at time t = we choose randomly the direction of time (increasing with probability 
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1/2) measuring with imaginary values. On each branch of the time axis we consider a Brownian 
motion, that is a process with independent increments over non-overlapping time intervals. To the 



(4.1) 



vector I -F(fi), . . . , F(t n ) j with < t± < . . . < tj < . . . < t n we associate the signed measure 
fi{F{ti) G dx u ...,F(t n ) G dx n } 

1 " 



Pr{B{ Sl ) G dx u ...,B{s n ) G dx n } + Pr{B{ Sl ) G dx u . . . , B(s n ) G dx n } 



n 



g 2i(t j -t j _ 1 ) 



'3— »*J 



3 tLj 



i y/2ni(tj -tj-i) fJi v /-27ri(t 3 - <j_i) 



da;. 



(27T) 



-n/2 



2 n,=i - *j— i 



= - z 4 n p ^3=1 2<(* i -t :f _ 1 ) , g*f»g -2«(*i-tj-l) 



J 3=1 



(27T) 



-n/2 



n™=i v*3 - *j-i 

We can verify that 



cos 



^ 2fo - 4 



3=1 



/ ... / n{F(ti) e d Xl ,...,F(t n ) e dx„} = l 



by writing (4.1) as 

(2tt)-"/ 2 



Vl 

E 



E 

vi=i 



(Xj - Xj-l) 2 , ,,7T . 

—f r — [n — 1) — I cos 



/ (X n - X n -l) 2 _ 7T 

\2{t n -t n ^) 4 



sm 



(a;„ - x ra _i) 2 _ 7r 
2(t„-t„-i) 4 



3=1 



(4.2) 



By integrating (4.2) with respect to x n the second term vanishes while the first one by applying the 
Fresnel integrals reduces to 



(27T) 



-(n-l)/2 



' n-l 



n ? = l 



COS 



(4.3) 



3=1 



The density (4.3) is the marginal Ai{n™ =1 -F(ij) € efej} of the vector 

(F{ti), . . . , F(t„_i)) . 
The one- and two-dimensional densities of the Fresnel pseudo-process read 

m/ \ i i dX\ ( x\ 7T 

fi{r (ti) G aiEi} = cos ' 

and 



T^rfT \2ii 4 



r / \ , , , , , dx\dxo ( x 2 (xn — X-i) 2 
H{F{H) G dxt,F(h) G dx 2 } = - .\ 2 . sin + ^ ^ 



(4.4) 
(4.5) 



By means of (4.1) we can construct the signed measure of cylinder sets of the form 



C = \ p| (a, < r, < h ; 
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in the following manner 



M {C} 



Ol 

bi 



/i{-F(*i) G <&!,..., F^) G dx n } 
(2tt)-"/ 2 



(4.6) 



COS 



(■Kj x j-l) 



E 



4 n d ^ 



The construction of (4.6) follows the same line of the signed measures of pseudo-processes related 



to higher-order heat equations (see for example, Krylov [15], Ladokhin [HHHI], Daletsky et al. 
Hochberg 10 ). The signed measure is extended to the er-algebra generated by the cylinder sets in the 
usual way (see, for example 16J ) . 



Remark 4.1. The signed measure constructed above is not markovian since 

n{F(h) G dx u F(t 3 ) G dx 3 \F{t 2 ) = x 2 } 

+ fi{F(h) G dxi|f (t 3 ) = a; 2 }/i{F(t 3 ) G dx 3 \F(t 2 ) = x 2 } 



(4.7) 



with ti < t 2 < t 3 . From (14. II) we have that 



MWi) G d Xl ,F(t 3 ) G dx 3 F(t 2 ) = ^2} = 



/x{F(ii) G dxi,F{t 2 ) = x 2 ,F(t 3 ) G dx 3 } 



n{F(t 3 ) G dx 3 } 



cos 



V 2|tl 



go) 1 (a 2 -a:i) . (x 3 -x 2 ) 
-tn| + 2|t 2 -ti| + 2|t 3 — t a | 



2tt 



rii=ife-*j-i) cos (! 



dx\ dx 3 . 



Clearly 



li{F(h) G dxi F(t 2 ) = x 2 }n{F{t 3 ) G dx 3 F(t 2 ) = x 2 } 



V2*i + 2|t 2 -ti| 2 J COS ^2t 2 + 



(X3 — X2) _ 7T 

2|t 3 -t 2 | 2 



V*l(*3-*2)(*2-*l) COS 2 (|| - f ) 



and thus we conclude that the statement (4.7) holds 



4.1 Feynman-Kac formula 

For a non-negative fc G C 2 (IR), the Feynman-Kac functional 

w(x,t) = E \ e -SS k ( F W)*>\F(0) = x 
must be understood in the sense that 

lim / ... / exp —/ > k(xj)(tj — ii-i) ^{-F(ii) G dxi, 



,F(t n ) G dx„} 



(4.8) 



(4.9) 



where \i is the signed measure defined in (4.1 ) and provided that the limit exists. This is similar to the 
definition of Feynman-Kac functional in Krylov [T^] , Lachal [T3] , Hochberg [TO] . We now show that 
the functional w solves the p.d.e. 



d 2 w 



(M) 



d 4 w 



d 2 w 



(4.10) 
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In order to prove (4.10) we consider that the measure [i appearing in (4.9) permits us to write w(x,t) 
as 



w(x, t) — — [w(x, it) + w(x, —it)] 



where 



w(x, t) = E 



-S*k{B{s))ds 



B(0) = x 



is the classical Feynman-Kac function for Brownian motion. Since 

dw 1 d 2 w 

= 2 'dx 2 V 33 '*) ~ Kx)w{x,t) 



by deriving (|4.10|) with respect to t we have that 

i d 2 i 
~2~dx 



OM) =o"5Z2"( a; ' it ) ~ ik(x)w(x,it) 



at 



i 2 w 

-^■(i, —it) + ik(x)w(x, -it). 



One more time-derivative in (4.12) yields 
d 2 w 



( )t^ t)= l 



d 2 fld 2 w 



dx 2 \ 2 dx 2 



[x, it) — k(x)w(x, it) I — ik(x) 



d 2 (Iff 



dx 2 \ 2 dx 2 



(x, —it) — k(x)w(x, —it) I — ik(x) 



1 d w 
2~dx^ 



2 dx- 



(4.11) 



(4.12) 



(x, it) — k{x)w{x, it) 



-{x, —it) — k(x)w(x, —it) 



We observe that for equation (4.10) a decoupling similar to that applied before in Section 3 and 4 does 
not work because of the form of the first-order time derivative (4.12). 



4.2 Superposition of vibrations 



We now analyse the superposition of solutions of the equation (2.1 ). Our approach is based on Fourier 
transforms and gives 



n 1 

n4 



+oo -i 



/2nt 



2t 4 



/ Xj 7T 
COS I —r~ - 

27rf \ 2t 4 



dxj 



dx 
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p 2 t 



( e 2 + e 2 J 
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2 n-l 
1 

Oll-l 



e i^-n _|_ e -j±^-n 



E 

fe=i 



n 
k 



1, 



e -iV(«-2) 1 



For n € 2N we obtain that 



cos 



n P 2 t_ 1 



i/2-l 



i E 

fc=0 



n— 2 / 
fc=2 v 



2 t(n-2k) 1 f n 



,-iV("-2t) 



2™ Vn/2 



(4.13) 
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whereas, for n £ 2N + 1 we have that 



1 ( \^ /2 fn\ (3 2 t(n - 2k) 

2^ D UJ C0S — 2 — ' 

k=Q v 7 



(4.14) 



Formula (4.13) shows that for n even the resulting superpositions of waves consists of components of 



the form (1.5) plus a Dirac delta function. 



Remark 4.2. For an even number of terms in the sum 53?= i ^jCO we nave delta components in zero. 
This surprising fact can be also confirmed by considering the convolution of two Fresnel waves as shown 
below. 
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4t V Vt 
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dw 



1 

4V7rt 
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it 



+00 



(u> 2 ) dw + 1 



At 
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V2 U V2 4t 
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2 V47rf V 4t 
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This result accords with (4.13) for n = 2 after an inversion of the Fourier transform. 



5 Vibrations of plates 

5.1 Vibrations of infinite plates 

The G?-dimensional version of equation (14. 71) has the form 



Ot 2 



1 

2 2 \dx 2 



(5.1) 
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and emerges in the study of vibrations of rigid thin structures like plates. The solution to (5.1 1 subject 
to the initial conditions 



u(xi,...,x d ,0) = l\ j=1 S(xj) 
Utixx, ...,x d ,0)=Q 



has the form 



i(xi, ...,x d ,t) 



cos 



J =1 



(5.2) 



(5.3) 



The equation of vibrating plates is examined in the book by Courant and Hilbert page 307] where 
the case of circular plates with Neumann boundary condition is outlined. 



Theorem 5.1. The Fourier transform of (5.3) is 



U(f3 1 ,...,f3 n ,t) = / e l ^=i^ x J U (x l7 ...,x d ,t)dx 1 ...dx d = cos ( V ) 



Proof. Let us write (5.4) as 
1 



jL sin - n ^ ./ 



e ipdXd cos 



j- 



2t 4 



(5.4) 



dx d 



■ ( x d A A 
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\2t 4 / 
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f ,<>:" ; i • 
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/2?rt 



+ 00 



x d_l 

2t 4 



— I dxd = cos '"" 



fit 



+00 



^-^U^sin^ 

2t A 2 
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H 



(5.5) 
(5.6) 



After m integrations, at the (m + l)-th step, we have that 
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^-^ d — rn — 1 2 
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2t 
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By performing the <i — m — 1 remaining integrals with respect to x±,.. . ,Xd- m -i we arrive at result 
(5l}. □ 
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In particular, for P± = . . . = f3 n = j3, from (5.4) we have that 

U (/?,..., p, t) = / e i ^t=^u(x 1 ,...,x d ,t)dxi...dx d = 



cos 



P 2 dt 



(5.7) 



which shows that the sum of the marginals of (5.3) has the same form of (1.5) at time td. For d = 1 



obtain the interesting result 
and the corresponding structure is depicted in figure [4j 



the equation (|5.1|) coincides with (|2.1|, (|5.3|) becomes (|1.5|) and (|5.4| reduces to (|2.9|). For d — 2 we 

1 



u(xi,x 2 ,t) = — sin 
Znt 



Xi ~\~ X<2 



2t 



(5.8) 




Figure 4: The vibrating surface (5.8) (within a sphere) at time t = 



It is apparent from (5.8) that the bivariate signed measure does not factorize. We can check that 



(5.8) (as well as (5.3)) integrates to unity 
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It is worthwhile to note also that 



u(xi,t) u{x2, t) u(xi, X2, t) 



because 



u(xi,t) u(x 2 , t) =— — COS 
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5.2 Vibrations of circular plates 

For the analysis of the vibrations of circular plates we need the following results from the theory of 
reflecting Brownian motion inside a circle Cr. Since the problem has isotropic structure (the circular 
invariance is due to the initial disturbance concentrated in the center starting off the vibrations) we 
can restrict ourselves to the Cauchy problem for the heat equation 



du 
~dt 



d 2 1 d 

dr 2 2 dr 



u, <r < R, i > 



(5.9) 



with boundary and initial conditions 



u(r,0) = S(r) 



du 



dC R 



= 0. 



(5.10) 



The solution of (5.9 1 with (5.10) is 



q ref (r,t) = 



(5.11) 



and is constructed by means of the inversion of radius. The probability law of the reflecting Brownian 
motion in Cr solves the adjoint equation 



du 
dt 



d 2 Id 1 

dr 2 r dr r 2 



(5.12) 



with initial and boundary conditions (5.10) and has the following explicit form 

r-2 R 4 R 4 



P re] >,<) = -e" 



i? > r > 0, t>0 



(5.13) 



and can be checked that it integrates to unity. The results above permit us to solve related boundary 
value problems for the equation of vibrations of plates. We have the following result. 

Theorem 5.2. The Cauchy problem 



dt 2 



1 



d 2 | 1 d 

dr 2 r dr 



u, < r < R, t > 



(5.14) 



with initial and boundary conditions 



has solution 



while the adjoint equation 

d 2 u 

W 2 



' u(r,0) = 5(r) 
u t {x,0)=0 
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t 2t t IrH 
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r f)2 



dr 2 



ld_ 
r dr 



-i 2 



u, <r < R, t> 



(5.15) 



(5.16) 



(5.17) 



with initial and boundary conditions (5.15) has solution 



r , V , T R . R 

p 1 (r, t) = - sm - — I — 5-sin— 5-. 
y V ' t It tH 2r 2 t 



(5.18) 
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Proof. We give the proof of (5.18) since the same method easily leads to (5.16|. The equation appearing 

u = (5.19) 



in (5.17) can be decoupled as 

d i ( d 2 Id 1 



d i ( d 2 Id 1 



dt 2 V dr 2 r dr 



and each Schrodinger-type equation is formally similar to (5.12) and we can thus write that 



P ref (r,t)=\ [p ref (r,it)+p re f(r,-it)] 
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□ 



The isotropic structure of the solution permits us to write the form of the circular vibrating plate 



as 



p(r,e,t) = —p ref (r,t) = — 
Fy ' ' ; 2n V ' 7 2nt 



r 2 R 4 R 4 
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2t r 3 2r 2 t 



< r < R, < 6 < 2w. 



(5.20) 



By considering the change of variable r' — R 2 jr in the following integral we get 
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Remark 5.1. The expression (5.20) in cartesian coordinates reads 

R 4 



p(x,y,t) 



1 

2~wt 



. x 2 + y 2 

sm 

2t 



R 4 



(x 2 + y 2 ) 4 2t(x 2 + y 2 ) 2 



(x,y) &C R ,t>0 



while the expression of q(r,9,t) = ^q re f(r,t) in cartesian coordinates reads 



q(x,y,t) 



2nt ^J x 2 + y i 



x 2 + y 2 
sin — — h sin 



R 4 



2t{x 2 +y 2 ) 2 



(x,y) €C R , t> 0. 



(5.21) 



(5.22) 



The cross-section of the vibrating disk (formula (5.16)) is depicted in figure [5] in three different 
instants. It shows that the vibrations initially extend to the whole plate and they fade off outside a 
neighbourhood of the origin as time passes. 



6 Subordination of the Fresnel pseudo-process with different 
processes 

The composition of the Fresnel process F(t), t > with different processes produces interesting results. 
The first one has been outlined in the previous sections and leads to the following statement concerning 
pseudo-processes of fourth-order which have been dealt by several authors and from different viewpoints 
(see Hochberg [TU], Nikitin and Orsingher [22], Nishioka [23], Lachal [HDl], Bcnacn °ur ct al. [3J). 
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Figure 5: The profile of the function (5.16) with < r < R = 1 and t = 1, 50, 100. 



Theorem 6.1. The pseudo-process F(\B(t)\), t > has measure density which satisfies the fourth- 
order heat equation 

u(x,0) = S(x). 1 ' 



Proof. We write 



= 2 



X 7T \ e 2t 



2tts \2s 4/ V2~7rt 



= ds. 



(6.2) 



The Fourier transform of (6.2) reads 



e if3x q{x,t)dx = I e 
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ds = cxp ( -— t ) . 



Clearly (6.3) satisfies the Fourier transform of the equation (6.1). Result (6.3) can also be obtained 
by exploting (2.9) and thus 



e^ x q(x,t)dx =2 
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If B(t), t > has volatility equal to a 1 then the composition F(\B(t)\), t > 0, has Fourier transform 
equal to 



ij3x 



q(x, t)dx = exp 



2 :! 



and for a =2 leads to the pseudo-process with density (6.2 1 



Remark 6.1. This result has been obtained in different ways in OUT]. Result (6.2) shows that the 
fundamental solution to the fourth-order equation (6.1 1 can be viewed as the profile of a vibrating rod 



with a Gaussian weight which damps the oscillations, see figure [6] 




Figure 6: The solution to the fourth-order heat equation (6.1) 



We examine here the composition of the Fresnel pseudo- process F(t), t > with the first passage 
time T t , t > (independent from F) of a Brownian motion B. The composition F(T t ), t > has the 
relevant property that its one-dimensional distributions are true probability distributions as the next 
Theorem shows. 



Theorem 6.2. The r.v. F(T t ) has density 

Pr{F{T t ) e dx}/dx = 



and (6.4) solves the fourth- order equation 
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satisfies equation (6.5). Since 
we have that 
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We now prove result (6.4) 
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Remark 6.2. We can confirm result (6.4) by means of a different approach. Since equation (6.5) can 
be written as 

d 2 .* <9 2 \ / <9 2 • * a 2 \ 

e-^)(^ + e--^)^0 (6.7) 



dt 2 



dx 2 



each component of (6.7) can be reduced to a Laplace equation by means of the time transformation 
t' = e *t. Therefore, the solution to (6.7) can be organized as 



Pr{F(T t ) G dx}/dx = 



1 
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(te^) 2 +x 2 {te-^Y+x 2 



t t 2 



nV2t 4 + x 4 



(6.8) 



This is tantamount to considering a Cauchy process C whose time either flows on the positive or 
negative imaginary axis. The direction of time is initially chosen with equal probability. 

Theorem 6.3. For the n-th order Fresnel iterated pseudo-process 

F n {t) = F 1 {\F 2 {\...F n+x {t)...\)\), t>0 
the governing equation of the measure density 



n{T n (t) G dx}/dx = 2' 



TT^±^cos \^-- A 



(6.9) 
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Figure 7: The density of the double-Cauchy r.v. F(T t ). 



with x = S\ and t — s„+i, satsfies the equation 

d 2 u ' - ^ 2 " +2 



(x,t) = ^2- 2 ^ +1 -^^^(x 7 t), x e E, t > 



subject to the initial conditions 



u(x, 0) = S(x), 
u t {x,0) = 0. 



(6.10) 
(6.11) 



Proof. One can prove the statement of the Theorem by evaluating the Fourier transform of (6.9 1 by 

yielc 

ds n 



successively applying formula (2.9). The Fourier integral w.r.t. x yields 
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The integral w.r.t. si yields 
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By iterating this procedure we arrive at the final relult 
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With this at hand result (6.10) immediately follows. 



□ 
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Figure 8: The table of fundamental solutions of the Fresnel fractional equation (3.1 1 with different 
orders 2v of fractionality. 
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